Anomalous diffusion in disordered media and random quantum spin chains 
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Using exact expressions for the persistence probability and for the leading eigenvalue of the Focker- 
Planck operator of a random walk in a random environment we establish a fundamental relation 
between the statistical properties of anomalous diffusion and the critical and off-critical behavior 
of random quantum spin chains. Many new exact results are obtained from this correspondence 
including the space and time correlations of surviving random walks and the distribution of the 
gaps of the corresponding Focker-Planck operator. In turn we derive analytically the dynamical 
exponent of the random transverse-field Ising spin chain in the Griffiths-McCoy region. 
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Ultraslow dynamics is a common feature of low- 
dimensional systems with quenched disorder in particular 
in the vicinity of a critical point. One of the well known 
examples in this respect is the one-dimensional diffusion 
process in a random media, when - in the absence of an 
average drift Srw - the mean-square displacement grows 
very slowly like jy 



\n 4 t 



(1) 



in contrast to the normal diffusive behavior (X 2 {t)) ~ t 
in the homogeneous case. The diffusion process remains 
anomalous for sufficiently small average drifts < Saw < 
when the average displacement has an algebraic 



time dependence: 



[<*(*)>]av ~ f 



(2) 



where the exponent < fi = [i(8rw) < 1 is a continuous 
function of the drift || . 

Another class of systems with ultraslow dynamical 
properties is represented by random quantum spin chains 
at very low temperatures. For example the asymptotic 
decay of the zero-temperature (imaginary time) autocor- 
relation function G(t) = [{erf {t)oj)] av at the quantum 
critical point (S — 0) is given by 



G(t,S = 0) ~ [lnt] 



-2x„ 



(3) 



where x m is the anomalous dimension of the average mag- 
netization. Away from the critical point, in the Griffiths- 
McCoy region with < S < 5 G the decay of the 
autocorrelations is of a power-low form 



G(t,S) -r 1 '* 



(4) 



where the dynamical exponent z(S) > 1 is a continuous 
function of the quantum control-parameter S. 

Comparing the basic dynamical properties of random 
walks in disordered environments and of random quan- 
tum spin chains one can easily notice close similarities, 



which hold both in the critical and off-critical situa- 
tions. A connection between the critical properties of di- 
rected walks and quantum spin chains is known for quite 
some time ||. It was also demonstrated recently that 
many surprising properties of the one-dimensional ran- 
dom transverse- field Ising model (RTIM), a prototype of 
random quantum spin chains, can be obtained very sim- 
ply through random walk arguments 0. 

In this Letter we go further and show that behind the 
similarities observed before there is a deep connection 
between the statistical properties of anomalous diffusion 
and the critical and off-critical behavior of the RTIM. We 
demonstrate this relation by comparing exact expressions 
for the random walk (RW) and that of the RTIM. In 
particular we show that the persistence probability of 
the RW and the surface magnetization of the RTIM have 
analogous forms and that the expressions for the leading 
eigenvalue of the Focker-Planck (FP) operator of the RW 
and the gap of the Hamiltonian of the RTIM are closely 
related to each other. We use then this correspondence 
to obtain new exact relations for the two systems, among 
others we present analytical results about the dynamical 
exponent z(S) in (|J). 

We start by considering the one-dimensional random 
walk with nearest neighbor hopping, which is character- 
ized by the transition probabilities it?i,i±i = w(i — ► i ± 1) 
for a random walker to jump from site i to site i ± 1. 
Here we are particularly interested in the general case, 
in which the transition probabilities are not necessar- 
ily symmetric ||, i.e. ^ Moreover, the 
random walker is confined to a finite number of sites 
i = 1, . . . , L. At the two ends of this interval, i.e. at 
i = and i = L + 1, we put adsorbing walls, which is 
simply modelled by setting wq_\ — wl+i.l = (i.e. the 
walker cannot jump back into the system once landed on 
or L + 1). The time evolution of the probability dis- 
tribution of the walk Pij(t), which is the probability for 
the walker to be at time t on site j once started at time 
on site i, is fully determined by the Master-equation 



1 



j- P.i(t) = M • Pj(t) ■ (5) where Z? = L J2i=i w i l+i 1S tne diffusion constant 

[ ]TT| | . Thus the average persistence probability in the ran- 
Here P t {t) = (P lfi (t), P hl (t), . . . , P ijL {t) , P %L+1 {t)) T and dom symmetric case scales as [p^™] ^ oc [D] av /L, sim- 
thc transition matrix or the Focker-Planck operator is ilarly to the homogeneous case. From here on we use 
{M)i,j = w iy j for i ^ j and (M)i,i = - J2j w i,o while the [• ■ .]«v to denote average over quenched disorder. 

In the general case, with non- symmetric transition 
probabilities we define the control parameter: 



initial condition is Pjj (0 



of the FP operator in (H) is defined by 



The eigenvalue problem 



UkK=u k 



(6) 



and all physical properties of the model can be expressed 
in terms of the left and right eigenvectors u k and v k , 
respectively, and the eigenvalues A&. For instance the 
probabilities Pi,j(t) are given by 



k 



Uk(i)vk(J) exp(Afet) 



(7) 



With adsorbing boundaries the two leading eigenvalues 
are zero and the corresponding eigenvectors are Vi(i) = 



ui(0) = 
ui(i) = ui(l) 



(8) 



i—l j 
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2,3,. 



while for the other zero mode V2(i) = V\(L + 1 — i) and 
similarly 112(1) =u%(L + l—i). The value of u±(l) in (||) 
is fixed by the normalization condition Ui(L + 1) = 1. 

We consider first a quantity that gained considerable 
interest recently in related models for anomalous diffu- 
sion PJlO[|: The persistence probability P pr (L,i), which 
is the probability that a walker starting at site i = 1 
does not cross its starting point until time t. Work- 
ing with adsorbing sites at i = and i = L + 1 
we have P pl (L,t) = Pii+i(t) and its long time limit 
Ppi(L) = lim^oo P pl (L,t) is simply given, via (0), 
Ppr(L) = ui(l)vi(L + 1), where we used the fact that 
there is no contribution from the second zero mode, since 
v 2 {L + I) = 0. Now with eq(||) we have the simple exact 
relation: 



Ppr(L)= 1 + EIT 



W 3,J-1 



i=l j=l J^ + L 



(9) 



Note that p pl (L) is the total fraction of walkers adsorbed 
by the right wall (i=L + 1) without ever having crossed 
the starting point. 

In a homogeneous medium with w^i+i = "U^+i.i = 
const we have p p ° m = (L + whereas in an inho- 

mogeneous environment with symmetric transition prob- 
abilities w iti+ \ = W l+1 .i 



[lni«<_] av 

[In io_»] 8 



var[lm 



var[ln"> 



(11) 



where (tc<_) stands for transition probabilities to the 
right (left), i.e. 10^+1 (w^i-x). For i5rw < there is an 
average drift of the walk towards the adsorbing site at 
i = L + l, therefore the persistence will have a finite value 
in the large system limit: lim^oo [p pr (L, S RW )] av > 0, 
whereas it goes to zero for <5rw > 0. 

Before we proceed with the analysis of the persis- 
tence probability (j^) we derive a similar formula for the 
largest non-zero eigenvalue of the FP-operator, the ab- 
solute value of which we denote by A m i n . According 
to the relation in (Q) the time-scale t r of the diffusion 
process is set by t r ~ A~[ n . It is technically easier to 
estimate A m j n using mixed boundary conditions, which 
will, however, not change the scaling behavior of A m i n : 
at i = we assume an adsorbing wall as before, whereas 
at i = L we impose a reflecting boundary by setting for- 
mally Wl,l+i = 0. Now, due to different symmetry of the 
problem there is only one zero mode of the FP-operator, 
and the second smallest eigenvalue in modulus will be 
Amin- To determine A m i n we use a perturbational method. 
First, we express the eigenvalue problem in eq(j(^) as 



- u(i)w it i-i + u(i + 1)10^+1 = —u(i)X T 



(12) 



where = and u(i) — u(i) — u(i — 1) ,i = 

1,2, ... , L in terms of the components of the left eigen- 
vector u = u min and u(0) = 0. Then, keeping in mind 
that we are interesting in situations when A TO j„(L) is a 
rapidly vanishing function of L we neglect the r.h.s. of 
eq([l2|) and derive an approximate expression for the left 
eigenvector from the first L — l equations of (|lj). Using 
this result we obtain an estimate for A m i n from the last 
equation of (|l2|): 



L-l 



" ,; 1 1 „.. 1 r "'./,, 1 



■ T CWL,L-1 = —7JT W L,L-1 ■ TT 

u(L) u{L) " Wj, j+ i 



, (13) 



which can be transformed into the final form by noticing 
that u(l)/u(L) = Pp r (L) in eq(|): 



L-l 



Amin ~ Pt>t{L)wl,L-1 ' TT '"' 

fJi W J.] + 1 



(14) 



L 



jj The scaling properties of A m in in (^4|) as well as the 

— , (10) persistence probability p pr (L) in (^) now can easily be 
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derived by establishing a correspondence of these quan- 
tities with the energy gap and surface magnetization, re- 
spectively, of the random transverse-field Ising model in 
one dimension defined by the Hamiltonian: 



H 



L-l 

E ■ l ' ,r '' T '- 1 

i=l 



L 

£ 

i=l 



(15) 



Here the af, cr| are Pauli matrices at site i and the Ji ex- 
change couplings and the hi transverse fields are random 
variables. The RTIM in ([l5]) has received much attention 
recently 12-14 |,[l5|,[7j and as was shown in there are 
simple expressions for the surface magnetization m s (L) 
as well as for the gap of the Hamiltonian A(L) of the 
RTIM in terms of the couplings and fields. Comparing 
those with our results in eqs(^) and (|lj) we can set up 
the following correspondencies 



(16) 



Consequently similar relations hold for the average quan- 
tities, when the transition probabilities (or equivalently 
the fields and the couplings) follow the same random 
modulation. In the following we use the correspondencies 
in (16) to derive new results. 

i) At the critical point with few — 0, which corre- 
sponds to the Sinai's walk jl[ , the distribution of the lead- 
ing eigenvalues of the FP-operator is very broad, A m j n (P) 
scales according to 



Wi,i+1 — 


-> J 2 


Wi,i-1 - 


h 2 


few — 


-» S 


Ppv(L) - 


- m 2 s (L) 


A m in(£) — 


- A(L) 


Mfew) " 





Amin(i) ~ exp(— const • L 1 / 2 ), 



TRW 



, (17) 



similarly to the analogous result for the energy gap of 
the RTIM 0. Note that this scaling relation |l|), which 
is consistent with the known relation between relevant 
time- and length-scales Q: L ~ (logi) 2 , can be most 
easily demonstrated by considering the probabilty distri- 
bution -Pl (In A m i n ), which is then expected to scale like 



Pi (In A min ) 



i5(lnA min /P 1 / 2 ) 



(18) 



as we confirmed numerically. 

ii) The scaling behavior of the persistence probabil- 
ity (for zero drift few = 0) in (||) follows also directly 
from the analogous result for the surface magnetization 
[m s (L)} av of the RTIM (?]]. Here we just have to men- 
tion that m s (L) at the critical point is not self-averaging, 
its average value is dominated by the rare events, which 
are of order O(l). From this follows that the same rare 
events determine the average of m 2 (L), thus the scaling 
behavior of [m s (L)] av and \m 2 (L)] av are identical. Then, 
using the correspondence in (|16|) we have the exact result: 



iii) In the non-critical situation with few ^ there is 
an average drift of the walk towards the site i — L+l (i = 
0) for <5rw < (few > 0). In the latter case, analogously 
to the surface magnetization of the RTIM , the average 
persistence probability vanishes exponentially for large 
system sizes 

[Ppr(i)]av ~ exp(-£/f), £ ~ few 

If the average drift of the walk is towards the adsorbing 
site at i = L + 1, thus few < 0, there is a non- vanishing 
infinite system size limit of the persistence probability 
(similarly to the existence of a finite average surface mag- 
netization of the RTIM @): 

(21) 



few > . (20) 



lim [p pr (L)]av ~ (-few)' 9pr , few < 



with /3 pr = 1, which is approached via an exponential 
size dependence. The corresponding correlation length 
is again given by £ ~ (—few) 2 , similar to the case 
few > in eq(j20|). Thus we can conclude that corre- 
lations defined on persistent walks are characterized by 
the average critical exponents 



1/2, 



2, /3pr 



(22) 



which satisfy the scaling relation /3 pr = 6v. 

iv) The time-dependent persistence 

probability P pr (L,t) introduced above is simply given by 
P pi (L,t) = Pi, L +x{t) = EfcWfc(l)w fc (P + l)exp(Afci). In 
the random, asymmetric case one expects the scaling re- 
lation 



[Ppr(LMt)l 



P pi (L/b,lnt/b^ 2 ) 



(23) 



when lengths are rescaled by a factor b > 1 and the re- 
lation in (|l|) or ( fT?! ) between time- and length-scales are 
used. Now taking 6 = L in the limit t — * oo we re- 
cover the exact result in (|l9|), on the other hand with 
b = (lni) 2 we have in the large system limit an ultraslow 
decay lim L ^ 00 [P pr (L, t)] av ~ (In*) -1 . 




[p pr (P)]av oc L- 9 , 8 = 1/2 few = 0. 



(19) 



0.5 1 1.5 2 2.5 3 3.5 

lnt/L 1/2 

FIG. 1. Scaling plot of the time-dependent survival prob- 
ability P pT (L,t) according to ( pj| ) for the asymmetric hopping 
model with a uniform distribution of hopping rates (averaged 
over 10 5 samples). The inset shows the corresponding scaling 
plot for the homogeneous case. 
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In the intermediate situation with kit ~ L 1 ^ 2 we have 
(with b = L) the finite size scaling form [P pi (L,t)] av ~ 
L _1//2 p(ln</i 1 / 2 ) with limy-Kx, p(y) = const — in con- 
trast to Pp° m (L,t) ~ L^^^t/L 2 ) in the homogeneous 
case. In Fig. 1 we show corresponding scaling plot for 
numerically generated results for finite systems that con- 
firm this scaling picture. 

v) Away from the critical point we reach the region of 
anomalous diffusion, which is equivalent to the Griffiths- 
McCoy phase of the RTIM. The relevant energy scale 
Amin(-k) for the RTIM) has a power law scaling 

behavior: 



Amin(i) 



-l/z(S) 



(24) 



and according to eq([l6|) the two exponents, // and l/z 
correspond to each other. At this point we use the re- 
sult that the value of [i is known exactly from the time 
dependence of the average desplacement of the walk in 
eq(^|) in the form B: 



1 . 



(25) 



Essentially this follows from the observation that for any 
independent identically distributed random variables x 
the distribution P(X) of A = X1X2X3 ■ ■ ■, which is remi- 
niscent of eq(|l4|), has an algebraic singularity at A = 
P(A) cx A~ 1+ ^ with fi given by [a^] av = 1, see @. 

Consequently we obtain for the dynamical exponent z 
of the RTIM in the Griffiths-McCoy phase the implicit 
equation 



l/z 



= 1 



(26) 



Note that for any distribution of J and h one obtains 
immediately the result l/z — 25 + 0(5 2 ), with 5 <C 1 as 
in ([n]), as has been observed earlier [ r2]|l5| ,|7j] . However, 
the exact result ( p6| ) is not restricted too <C 1, but is 
valid in the whole Griffiths-McCoy region. For example 
for the uniform distribution 

tt(J) = 6(1 - J)e(J); p(h) = h^e(h - h)Q(h) , 

the dynamical exponent is given by the solution of the 
equation 



zlog(l 



2 ) = -lnfr (=-25) 



(27) 



The relation (g7|) is indeed satisfied by the numerical es- 
timates for z reported in [p"3 15 7[|. 



To summarize in this letter we have revealed a funda- 
mental relation between the anomalous diffusion of ran- 
dom walks in disordered environments and the slow dy- 
namics, at criticality and in the Griffiths-McCoy region 
of the random transverse Ising chain. With this analogy 
at hand we were able to derive a number of new exact 
results for both systems. Many new applications of the 



above mentioned analogy are obvious: there is an enor- 
mous number of exact results for various quantities of 
random walks in random one-dimensional environments, 
and most probably many of them can be directly trans- 
ferred to corresponding quantities of random quantum 
spin chains near the quantum critical point. It remains 
a subject of future research to study in how far these 
relations carry over to higher dimensions. 
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